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Abstract 

Based on the Holstein double-exchange model and a highly efficient single cluster Monte Carlo approach we study 
the interplay of double-exchange and polaron effects in doped colossal magneto- resistance (CMR) manganites. The 
CMR transition is shown to be appreciably influenced by lattice polaron formation. 
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Over the last decade the magnetic and transport 
properties of CMR manganites (Lai_x[Ca, SrJxMnOa 
with 0.2 ^ X < 0.5) have attracted a considerable 
amount of research activity, and in particular pola- 
ronic features near the transition from the ferromag- 
netic to the paramagnetic phase remain to be an in- 
tensely studied subject. A realistic description of the 
observed Tc and of the electrical resistivity data is com- 
plicated by the requirement of incorporating strong 
electron-phonon interactions in addition to the mag- 
netic double-exchange (DE) [1] . In the present study we 
aim at developing realistic effective models for the spin 
lattice interaction and appropriate Monte Carlo (MC) 
techniques for their simulation. Our starting point is 
the Holstein-DE model 

H = -Yl *'J 4 ^j ~ V'^p'^o ^{bl+bi) c\ c, +0^0 ^ h]bi , 

{i,i) i i 

(1) 

where the first term describes the well known DE 
interaction, characterised by the transfer amplitude 
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which 
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depends on the classical spin variables {9i,(j)i}. The 
second term accounts for a local coupling (oc Sp) of 
doped carriers to a dispersionless optical phonon mode 
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with frequency ujo, and the last term refers to the 
dynamics of the harmonic lattice. 

As a first step let us focus on the numerical solu- 
tion of the DE part, which is characterised by non- 
interacting fermions coupled to classical spin degrees 
of freedom. In a MC simulation of such types of models 
the calculation of the fermionic energy contribution, 
which depends on the classical degrees of freedom, is 
usually the most time consuming part, and an efficient 
MC algorithm should therefore evaluate the fermionic 
trace as fast and as seldom as possible. The first re- 
quirement can be matched by using Chebyshev expan- 
sion and kernel polynomial methods [2] , but so far this 
approach was combined only with standard Metropo- 
lis single-spin updates [3] . For the second requirement 
a hybrid MC involving classical time evolution of an 
effective spin model and an approximate diagonalisa- 
tion of the fermionic problem was suggested [4] . How- 
ever, both approaches have a few drawbacks, the flrst 
suffers from the frequent evaluation of the fermionic 
trace, the latter is rather complicated since it involves 
a molecular dynamics type simulation of the classical 
degrees of freedom. 

As an alternative we propose a combination of ker- 
nel polynomial expansion and cluster MC (CMC) [5], 
which is known to be the fastest approach to non- 
frustrated classical models. Averaging the first part 
of (1) over the fermionic hopping we can define an ap- 
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Fig. 1. Comparison of magnetisation data for the pure DE 
model at a; — 0.5, obtained witli the effective model (2), the 
hybrid MC approach [4], and the new CMC method. Here and 
in the following all energies are measured in units of t. 

proximate classical spin Hamiltonian 

HcB = "-^off ^ v 1 + Si ■ Sj , (2) 

which can be easily simulated with a rejection-free 
Wolff [6] single cluster type approach. Setting Jos — 
x{l — a;)/%/2, where x is the hole concentration, the 
magnetisation data and critical temperatures of this 
approximate model agree surprisingly well with the 
full DE system (cf. Fig. 1). Thus, the latter can effi- 
ciently be simulated by the new hybrid approach: (i) 
all spins of the system are altered by multiple clus- 
ter flips, where preferably Jett is random with (Jeff) = 
x{l~ a;)/\/2), (ii) the fermionic energy is evaluated us- 
ing the kernel polynomial method, and (iii) the new 
spin configuration is accepted with probability P{a — > 

b) = min (l, ZlTAtZt )- Here W{a), W{b) are the 
Boltzmann weights of the spin configurations a, b, and 
A[a -^ b), A{b —> a) are the accumulated a priori prob- 
abilities of the cluster flips connecting a and b [5] . Of 
course, the MC update is no longer rejection-free, but 
since Jfcff yields very good estimates for H, P{a -^ b) 
is sufficiently high for a reliable and fast simulation of 
reasonably large clusters. This allows for a confirma- 
tion and possible future extension of previous results 
for the DE model [4]. 

Encouraged by the above findings we can now rely 
entirely on the effective spin model Hcb and include lat- 
tice polaron effects on a variational level. Applying the 
modified variational Lang-Firsov transformation [7], 
from (1) and (2) we obtain the effective polaron-DE 
model, 
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+ JI^,+Ne,x[{l^^f{l~x)-l], (3) 

where 7 is a variational parameter which measures the 
importance of the polaron effect (0 < 7 < 1). The 



Fig. 2. Left: Magnetisation data (filled symbols) and polaron 
parameter 7 (open symbols) for typical values of Cp and uq. 
Right: Dependence of the critical temperature on the clcc- 
tron-phonon interaction strength. 

above model is treated by single cluster MC, where the 
optimal 7 is adjusted after each cluster flip. Note that 
the spin-lattice part of (3) is symmetric with respect 
to a; = 0.5, i.e., to explain the asymmetry of the man- 
ganite phase diagram additional Jahn- Teller type lat- 
tice interactions need to be included. Figure 2 shows 
the magnetisation as well as 7 for a few typical pa- 
rameter sets. Clearly the polaronic effect (larger 7) is 
most pronounced near and above the critical tempera- 
ture. As expected, the electron-phonon interaction re- 
duces the critical temperature, and for small phonon 
frequencies (adiabatic case) the ferromagnetic phase 
may cease to exist completely. Potential future exten- 
sions of the present work could include Jahn- Teller 
modes [8] and studies of the conductivity, e.g., along 
the lines of Ref. [9]. In addition, the formation of mag- 
netic polarons in the DE model should be accessible 
within the proposed CMC approach. 
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